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Abstract 

A  class  of  singular  integro-differential  equations  in  Lebesgue  spaces  are  stud¬ 
ied.  There  are  many  applications  of  the  singular  integro-differential  equations 
discussed  in  this  paper.  An  example  in  modeling  the  stress  distribution  of  an  elas¬ 
tic  medium  with  holes  is  discussed  in  the  paper.  Direct  numerical  schemes  using 
a  collocation  method  and  a  mechanical  quadrature  rule  designed  for  the  singular 
integro-differential  equations  are  proposed  for  arbitrary  smooth  closed  contours. 
Convergence  analysis  of  these  methods  are  given.  Numerical  examples  are  also 
provided. 

Keywords:  singular  integro-differential  equations,  elasticity,  collocation  methods, 
mechanical  quadrature  rule. 


1  Introduction 


In  this  paper,  we  study  the  following  system  of  singular  integro-differential  equations 
(SIDE) 


( Mx  =)  y ^[Ar(t)x^r\t)  +  Br(t) —  [  - - — dr 

7^0  mJ  r-t 

+  J  Kr(t,  t )  •  x(r)(r)dr]  =  /(f),  t  €  T, 


(1) 


where  Ar(t),  Br(t),  and  Kr(t,  r),  (r  =  0,  u)  are  given  m  by  m  matrix  functions,  f(t)  is  a 

given  m  dimensional  vector  function,  x^°\t)  =  x{t)  is  a  m  dimensional  unknown  vector 

dr  x(t^)  _ 

function  with  up  to  z/s  continuous  derivatives  —  (r  =  1  ,v),  and  T  is  an  arbitrary 
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smooth  closed  contour.  We  wish  to  find  the  solution  x(t)  to  equation  (1)  satisfying  the 
condition 

— —  [  x(T)T~k~1dT  =  0,  k  =  0,  v  —  1.  (2) 

2m  J 
r 

We  will  develop  two  direct  methods  to  solve  the  system  of  singular  integral  equations 
(SIDE)  (l)-(2)  and  provide  the  convergence  analysis. 

Singular  integral  equations  (SIE)  and  SIDE  have  been  used  to  model  many  physical 
problems,  for  examples,  elasticity  theory  [2,  9,  11,  15,  16,  17,  20],  aerodynamics  [11,  12, 
14] .  We  present  an  example  from  elasticity  theory  in  Section  2  to  show  one  application 
of  the  SIDE  (l)-(2)  discussed  in  this  paper. 

Since  analytic  solution  to  SIDE  is  rarely  available,  we  look  for  an  approximate 
solution  to  the  SIDE  using  direct  methods.  Most  of  early  direct  methods  for  SIE  and 
SIDE  are  designed  for  the  case  where  the  boundary  T  is  a  unit  circle  or  an  interval 
[1],  [3]-[7],  [8,  18].  The  requirement  of  unit  circle  is  essential  because  some  special 
polynomial  interpolations  are  only  defined  on  the  unit  circle,  see  [13],  for  example. 
There  are  few  papers  that  deal  with  arbitrary  closed  contours  in  the  literature. 

It  should  be  mentioned  that  using  a  conforming  mapping,  one  can  transform  an 
arbitrary  smooth  closed  contour  to  the  unit  circle.  However,  this  approach  may  not 
simplify  the  problem  due  to  the  following: 

•  It  is  not  an  easy  job  to  find  a  conforming  mapping. 

•  The  coefficients,  the  kernels,  and  the  right  hand  side  of  the  transformed  equation 
may  be  more  complicated  and  may  not  be  smooth  anymore. 

•  The  convergence  analysis  may  be  more  complicated  due  to  the  transformation  of 
the  contour. 

In  the  1980-90’s,  theoretical  foundations  for  direct  methods  for  SIE  and  systems 
of  SIE  defined  on  arbitrary  smooth  contours  in  the  complex  plane  were  established 
in  [22,  23],  and  maybe  others.  But  direct  methods  for  SIDE  or  system  of  SIDE  for 
arbitrary  contours  are  yet  to  be  developed  and  analyzed. 

In  this  paper,  we  develop  a  collocation  method  for  the  systems  of  SIDE  with  arbi¬ 
trary  smooth  contours,  and  the  mechanical  quadrature  rule,  see[23]  for  the  definition,  to 
approximate  the  integrals  involved  in  the  system  of  SIDE.  We  also  give  the  convergence 
analysis  of  our  new  method  in  Lebesgue  spaces. 

The  classical  continuous  function  space  can  not  be  used  because  the  singular  oper¬ 
ator  of  the  integration  is  unbounded.  The  classical  theory  of  projection  methods  for  Lp 
spaces  does  not  apply  because  the  norm  of  some  projectors  (for  example,  the  projectors 
of  interpolation)  is  unbounded  for  Lebesgue  spaces.  Thus,  it  is  necessary  to  elaborate 
a  new  theory  for  the  collocation  and  quadrature  method  discussed  in  this  paper. 
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The  paper  is  organized  as  follows.  In  the  next  section,  we  present  an  application  of 
SIDE  to  show  the  importance  of  SIDE.  In  Section  3,  we  describe  the  numerical  schemes 
for  systems  of  SIDE  with  an  arbitrary  smooth  contour.  The  convergence  analysis  is 
given  in  Section  4.  Numerical  examples  are  provided  in  Section  5. 


2  An  application  of  singular  integro-differential  equations 


We  show  an  application  of  SIDE  in  this  section  to  show  the  importance  of  our  work  on 
SIDE. 


Consider  an  elastic  medium  under  conditions  of  plane  strain  or  generalized  plane 
stress  which  fills  all  of  the  infinite  plane  of  the  variable  z  and  is  weakened  by  a  hole  of 
arbitrary  shape.  By  convention,  this  problem  refers  to  the  determination  of  the  stress 
components  along  the  contour  of  the  hole  when  external  forces  are  prescribed  on  the 
same  contour.  The  state  of  stress  at  infinity  is  given. 


As  described  in  [9],  the  physical  model  is  the  following: 

f  41)  +  ^D  =  2($l(a+^(0), 

|  <41}  -  cr(p]  +  2  MV  =  -L—  (w(£)$i(£)  +  w'(0^i(0)  > 

l  p2u  (£)  v  7 

where  cr^\  rW  are  the  components  of  the  stress  in  the  curvilinear  coordinates 
related  to  the  conformal  mapping 


Si(£)  = 


z  =  w(£),  £  =  pev 

Mo 


MO  = 


MO 


<«/(£)’  w  MO 
MO  =  </>M£)) ,  MO  =  'i/’MO)> 

and  <j>(z),  ip(z),  as  always,  denote  the  complex  potentials  in  the  z  plane. 


(4) 

(5) 


Suppose  that  the  region  E  with  boundary  7  which  is  mapped  by  relation  (4)  onto 
the  physical  region  S  represents  an  infinite  region  in  the  plane  of  the  auxiliary  variable 
£  outside  the  unit  circle  centered  at  the  origin.  This  means  that  the  function  cu(£)  is 
representable  in  the  region  E  as  the  series 


w(£)  =  J2£  +  j  +  p  +  .. 


e 


(6) 


It  may  also  be  recalled  that  the  functions  <^>i(£),  V’i(£)  admit,  in  the  same  region 
the  representation, 


MO  =  -0Xn?\  log£  +  RTZ  +  MO 

27t(1  +  K) 

MO  =  ~  t~t  log£ +  RT't  +  MO 

Z7T{l  +  AvJ 
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where  ip*  are  homomorphic  functions,  regular  at  £  =  oo,  and  X,  Y  are  the  com¬ 
ponents  of  the  resultant  of  the  external  forces  applied  to  the  boundary  of  the  region 
5. 

Assume  now 


<T(i)-<7+(To  (i)  ,  o  T(i)_T  ,To 

Up  —  U  p  ~r  Up,  Uv  —  Uy  T  Uy,  I  rv  —  !rv  T  lrv, 


(8) 


r°  a0  t° 

rv 

constants  T  and  T  (T  =  f),  and  represent  relations  (3)  with  p  =  1  as  follows  (a  =  ew ): 


where  07,  a®,  t®v  are  the  components  of  stress  in  a  uniform  field  characterized  by  the 


2  ($(cr)  +  $(<r))  =  (Tp  +  av  on  7, 

(w(ct)4>  (ct)  +  cj  (cr)'I'(cr))  =  cr2uj'(a)  (av  -  ap  +  2iTpv) 


(9) 


Here  $  and  'h  are  two  unknown  functions,  homomorphic  in  the  region  E.  From  (5)  to 
(7),  we  have  the  following  asymptotic  when  |£|  is  large, 


71  ...  A 


m  =  'i+o{c'z),  no  =  j+o(c2), 


A  =  - 


X  +  iY 
2irR(l  +  k) 


A'=  <X-iY ) 


2irR(l  +  k) 


The  stresses  a®,  a®,  t®v  then  can  be  calculated  from 


o°p  +  o°  =  2^0(0  +  M0), 


Ov-°p  +  2iTPv  = 


2£2 


(W(e)$o(e)+w'(o®o(o 


with 


4>n  = 


P2"'(0 

RT  RT' 

u  (0  u  (0 


Compared  with  (8),  we  have 

MO  =  M)  +  MO,  MO  =  *(£)  +  MO- 


(10) 

(ii) 


(12) 


(13) 


It  is  well  known  that,  whenever  the  region  occupied  by  the  elastic  medium  contains 
the  infinitely  remote  point  of  the  plane,  the  constants  T,  r\  X,  Y  are  specified  by 
the  problem  itself.  The  stresses  <7°,  a t®v  and  the  constants  A  and  A'  are  therefore 
known,  they  are  defined  by  formulas  (11)  and  (12). 

Our  problem  is  to  determine  the  normal  stress  av  by  formulas  (9)  from  the  given 
values  on  circumference  7  of  the  other  two  components,  ap  and  rpv.  Assume 

(Tv  +  (Tp  =  fi(o-),  — 2(crp  -  i.Tpy)  =  O (a)  on  7  (14) 
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and  rewrite  (9)  as 


2  ($0)  +  $(<r))  =  n{a), 

(t2lu' (a)  (Q(a)  +  Q(cr))  —  2lu((j)& (a)  =  2lu' (a)'k(cr)  on  7. 


(15) 


After  the  transformations  we  obtain  the  required  equation  for  determining  H(<r)  as 
follows: 


<j2lu' (a)[£l(cr)  +  0(<r)]  +  — 7  f  t2o/(t)[H(t)  +  0(r)] 

m  Jy 

1  /'  uj(t)  —  to  (a)  1  4 AR 

- /  - —  n  (r)dr - =  0  on  7, 

m  Jj  t  —  a  a 


dr 


t  —  a 


(16) 


where  A,  R,  0(<r)  are  known.  After  some  manipulation,  see  [10]  as  well,  we  obtain 


-y77Tn/  ,  cr2W  (<t)  /  n(r) 

<j2c 0  (<t)H(<7)  H - /  - -dr 

m  Jy  T  —  (7 

1  f  w(r)  -  w(a)n,.  ,  1  f  rV(r)  -  crV(cr) 


7TI 


7 


t  —  a 


-fi  (r)dr  H - / 

m  Jy 


t  —  a 


Q(r)dr  (17) 


4AR 


<7 


,  1  f  rV(r)0(r) 

—  cr-w  (<7)0(<j) - ;  /  - dr 


m 


t  —  a 


or 


A(fj)fl(fj)  +  dr  +  —  [  Kv (<7,  t)H^(t) dr  =  f  (a) , 

m  J.y  t  —  a  m  47 


where 


4/  4  2  7  1  7..'/  \  rV(r)  -  o-2w'(o-)  u(<t)-u(t) 

A{a)  =  a2uj  (a),  K0{<j,t)  =  - ,  Ai(cr,  r)  = 


t  —  a 


t  —  (J 


fi?) 


4AR 

a 


a2  u  (a)  0(a) 


jL  f  T2Uj\t)Q(t) 
ni  Jy  t  —  cr 


If  uj(a),  f(a),  A(a)  are  Holder  functions  on  7,  then  Ku(a,r)  =  ^J^x  ,  and  Ku(a,r), 
are  also  Holder  functions  for  0  <  A  <  1,  and  ^  =  0,1,  see  [10].  Thus  the  application 
problem  can  be  described  by  the  SIDE  (l)-(2)  with  A  =  0. 


3  A  collocation  and  quadrature  method  for  system  of 
SIDE 


Before  we  describe  our  numerical  schemes,  we  need  some  theoretical  preparations. 
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3.1  Theoretical  preparations 

Let  r  be  a  smooth  Jordan  curve  which  divides  the  entire  complex  plane  C  into  two 
parts,  F~  and  F+ .  Assume  that  F~  contains  the  infinity  and  F+  contains  the  origin. 
Therefore  F~  =  C  —  {F+  UT}.  Let  z  =  tp(w)  be  an  analytic  function  that  maps 
conformally  the  exterior  of  the  unit  circle  to  F~  so  that  'ip(oo)  =  oo,  V’ '(oo)  >  0. 

Let  us  consider  the  Riemann  function  t  =  ip(w),  w  G  To  =  {£  :  |£|  =  1},  that  has  at 
least  up  to  second  derivatives.  Assume  that  tp^2\w)  satisfies  the  Holder  condition  on 
|ir|  =  1,  that  is  ip(2\w)  G  H^Tq),  for  some  parameter  //.  p  (0  <  p  <  1).  We  denote 
the  set  of  such  contours  by  C( 2-,p),  see  [23]  for  more  details. 

In  the  complex  [Lp(T)]m(l  <  p  <  oo)  space,  the  norm  of  a  vector  functions  g(t)  = 
■  ■  ■ ,gm(t ))  is  defined  as, 


m 


llsll  =  E  II® Hr! 

k= 1 


llfffcllp  — 


where  l  is  the  length  of  T. 

We  will  develop  a  direct  method  for  the  systems  of  SIDE  (1)  together  with  conditions 

(2). 

3.2  A  collocation  and  quadrature  method 

Using  the  operators  P  =  ^(I  +  S),  Q  =  I  —  P,  where  I  is  the  identity  operator,  and  S 
is  a  singular  Cauchy  nucleus,  equation  (1)  can  be  written  as  follows: 

V 

(Mx  =)  Y  ^Ar(t)(Px^)(t)  +  Br(t)(Qx^)(t)+ 

r= 0 

+  2^i  J  Kr(t,T)x{r\'r)drSj  =  f(t),  ter  (18) 

where  Ar(t)  =  Ar(t)  +  Br(t ),  Br{t)  =  Ar{t)  —  Br(t),  r  =  0,  v. 

We  seek  an  approximate  solution  to  problem  (l)-(2)  of  the  form 

n  —1 

Xn(t)  =  Y  +  E  &)tk’  t  G  r’  (19) 

k= 0  k=—n 

where  =  £&>  k  =  —  n,  n  are  unknown  vectors  of  dimension  m.  Note  that  xn(t)' s 
constructed  using  (19)  satisfy  the  condition  (2). 
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In  our  collocation  method,  we  choose  2 n  +  1  different  points  on  T  {tj}  such  that 
the  SIDE  can  be  satisfied  exactly 


(■ Mxn){tj )  -  f{tj)  =  0, 


(20) 


As  a  result,  we  obtain  a  system  of  linear  algebraic  equations  for  the  coefficients 
k  =  —  n,  n, 

s  ir)rL w- 


+  S(  11  (*-D! 


1  _(fc  +  u) 

(h 

k= 0 

(A:  +  r  —  1)!  1 


■‘A'  '  «-*  +  2S  '  t  (iT^y  /  .6  (21) 


J  Kr(tj,T)T  k  rdT  ■  £_fc}  = 


where  j  =  0,  2 n. 


We  need  to  approximate  the  integrals  in  the  system  of  linear  algebraic  equations 
above  as  well.  This  is  done  using  the  mechanical  quadrature  formulae  described  in  [23] , 

hi  J«^T‘+kdT  ahS  [/”<T'+1  '  s(r)V-'dr. 

r  r 

where  the  operator  of  interpolation  Un  is  determined  by,  see  [23] , 


2  n 


W)  =  n 


( Ung)(t )  =  ^2g(ts)  ■  ls(t) 

s= 0 

2n  t-tk  (t'n 


k=0  ,k^j  tj  tk  \  t 


(22) 


t)  =  E  Ak)tk’  ter’  J  =  0, 2n. 


k=—n 


Thus,  the  system  of  linear  algebraic  equations  for  the  unknown  of  (21)  has  been 
approximated  by  the  following  system  of  linear  algebraic  equations 

UMtj)  E  „( 


r=0 


(k  +  v  -  r)!  J 


k= 0 


+ Mb)  E(-i  •  £-* 


fc=i 


(fc-1)! 


(23) 


/c=0 


5—0 


2n 


+  E(-ir(fc(^  1)])!Eifr(ii,t,)i.-1A'li->}  =  /(%), 


k=  1 


for  j  =  0, 2 n.  Thus  we  have  obtained  the  system  of  linear  algebraic  equations  for  the 
unknowns  £*.  in  (19).  Once  these  coefficients  are  obtained,  the  solution  to  the  singular 
integro-differential  equations  is  expressed  as  (19). 
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4  Error  estimates  for  the  numerical  method 


In  this  section,  we  give  some  error  estimates  for  the  numerical  method  proposed  in  the 
previous  section. 


As  was  proved  in  [10],  a  paired  vector  functions 


dT(Px)(t )  ^  du(Qx)(t) 


dty  dtv 

represented  by  the  integrals  of  Cauchy  type  with  the  same  density  ((t): 


can  be 


dv{Px){t) 

dty 

du(Qx)(t) 

dP 


T  —  t 


dr, 


C(t) 

T  —  t 


dr , 


t  G  F+ 
t  G  F~.  | 


(24) 


By  means  of  the  above  formulae,  the  problem  (l)-(2)  can  be  reduced  to  an  equivalent 
(in  terms  of  solving)  singular  integral  equations  system 

(RC  =)C(t)C(t)  +  L>(^  I  ^-dr  +  [  h(t,T)C{r)dT  =  f(t),  teT,  (25) 

m  J  t  —  t  Zm  J 
r  r 


where 

C[t)  =  -\Au(t)  + 1  vBu{t)}, 


(26) 


and  h(t,r)  is  a  matrix  function  belonging  to  [C(r)]m  space  for  both  variables  t  and  r, 
see  [10].  The  equivalence  of  the  existence  of  the  solution  between  the  system  of  SIE 
(25)  and  the  problem  (l)-(2)  is  the  result  of  the  following  lemma  from  [10]. 


Lemma  1  The  system  of  SIE  (25)  and  the  problem  (l)-(2)  are  equivalent  in  terms  of 
the  solvability.  That  is,  for  each  solution  £(i)  of  system  (25),  there  is  a  solution  x(t) 
of  the  problem  (l)-(2),  determined  by  the  following  formulae 

{Pxm  =  2XZ1) !  /  C<r)  (<T  “  1  l0g  (‘  -  ;)  +  E  “T  (2T) 

jp  \  k — 1  / 

(W(‘)  =  J'Td,  /  C(r)T-  ((r  -  t)-‘  log  (l  -  ()  +  E  Ar-*-v)  dr, 

where  a k  and  (3^,  are  real  numbers.  On  the  other  hand,  for  each  solution  x(t)  of  problem 
(l)-(2),  there  is  a  solution  ((t) 

_  du(Px)(t)  du(Qx)(t) 

U  ’  dtv  +  ’  dtv  ’ 

to  the  system  (25).  Furthermore,  given  a  set  of  linear-independent  solutions  £(i)  to  the 
system  (25),  there  are  corresponding  set  of  linear-independent  solutions  to  the  problem 
(l)-(2)  from  (27)  and  vise  versa. 


The  proof  can  be  found  in  [10]. 

Before  we  discuss  the  main  convergence  theorem  of  our  method,  we  need  several 
additional  lemmas  from  [19].  Those  lemmas  will  be  used  in  our  discussions  in  this 
section. 

Define  the  space 

[Wp)m(r)  =  {g-  3gW  G  [C(T)]m,  r  =  0^1,  gM  G  [Lp(T)]m}  . 

o  (y)  o  W 

For  any  g  G  \Wp  ]m,  g  satisfies  equation  (2).  The  norm  in  [Wp  }m  is  defined  by 

Ml  P,u  =  \\gM\\[Lp]m- 

We  shall  denote  the  image  of  space  [Lp]m  under  the  mapping  P  +  t~uQ  by  [LPt u\m  with 
the  same  norm  as  in  [Lp\m. 

o  (v)  .  . 

Lemma  2  The  differential  operator  Du  :  [Wp  ]m  — >  [Lp,u\m,  ( Dug)(t )  =  g^'(t)  is 

_  o  (D 

continuously  reversible  and  its  reverse  operator  D  v  :  [LPiV\m  — ►  [Wp  ]m  is  determined 
by  the  equality 

(D-“am  =  (iV+9)(t)  +  (JV-„)((), 

where 

(JV+9)(t)  =  JC-  1)!  /(Ps,(t)(t  “  *)^1  log  (‘  -  l)  dT • 

r 

{N~Sm  =  1)!  /(<39)(t)(t  “  l0g  (‘  -  9 

r 

The  proof  can  be  found  in  [19]. 

o  (v) 

Lemma  3  The  operator  B  :  [Wp  ]m  —■ y  [Lp\m,  of  the  form  B  =  (P  +  tvQ)Dv  is 
reversible  and 

B-1  =  D-v(p  +  t-VQy 
The  proof  can  also  be  found  in  [19] . 

The  existence  of  the  solution  to  the  systems  of  linear  algebraic  equations  (21)  and 
(23)  is  given  in  the  following  theorems. 

Theorem  4  Assume  the  following  conditions  are  satisfied: 

1.  the  contour  T  £  C{ 2,p),  0  <  p,  <  1, 
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2.  the  matrix  functions  Ar(t)  and  Br(t )  belong  to  the  space  [Ha(T)]m,  0  <  a  <  1, 
for  r  =  0,  v, 

3.  det(Au(t))  ■  det(Bu(t))  /  0,  for  any  t  G  T, 

4 ■  the  left  partial  indexes  of  tuBf1(t)Au(t)  are  zero,  for  any  t  £  T, 

5.  the  matrix  functions  Kr(t,r )  £  Hp[(T)]m,  0  <  (5  <  1  with  variables  t  and  t,  and 
the  vector  function  f(t)  £  [C(r)]m  for  any  t  £  T, 

o  M 

6.  the  operator  M  :  [Wp  }m  — >  [Lp(T)]m  is  linearly  reversible, 

7.  the  points  tj  (j  =  0,  2 n)  form  a  system  of  Fey er  knots,  see  [21]  for  the  definition, 
on  T: 


tj  =  il)(ex P  (^T[(J  “  n)))  ’  J  =  °>  2n> 


=  -1. 


(28) 


Then,  if  n  is  large  enough,  the  system  of  linear  equations  (21)  has  a  unique  solution  (k, 

k  =  —n,  n.  The  approximate  solutions  xn(t),  constructed  from  formula  (19),  converges 

o  (y) 

as  n  — ►  oo  in  the  normed  space  [Wp  ]m  to  the  exact  solution  x(t)  of  SIDE  (l)-(2). 
Furthermore,  the  following  error  estimate  holds: 


\x  —  x 


n  |  \p,v  O 


1 


nu 


+  O  (u(f-  * ))  +  O  (u\h-,  * ))  r=7  Sn.  (29) 


Proof:  We  show  that  if  n  is  large  enough,  then  the  operator  UnMUn  is  reversible. 
The  operator  acts  on  the  subspace 


f  n  —1 

[Xn]m=\t^aktk+  £  <*ktk) 


k= 0 


k=—n 


° 

with  the  norm  defined  in  [Wp  \m  to  the  subspace 

[Rn}m  =  |  E  rktk  | ,  ter, 

I k=—n  I 


with  the  norm  defined  in  [Lp(T)]m. 

Similarly  to  (24)  we  represent  ^  ^  ^  and  — by  integrals  of  Cauchy 

type  with  the  same  density  £n(t)  : 


dtv 


dth 


dv{Pxn)(t)  =  J_  r  Cn(r)dT)  t  €  F+ 


dtv 


2m  J  t  —  t 


du(Qxn)(t)_  =  t_^_  r  Cn(r)  t  ^ 


dth 


2iri  J  t  —  t 


(30) 
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It  is  easy  to  prove  that 


Cn(t)  = 

k= 0 


{k +  k 


k\ 


k= 1 


{k  +  u-l)\^_k 


(k-  1)! 


and  therefore  ( n(t)  £  [ Rn]m ,  for  all  t  E  T. 

Using  (30),  we  reduce  the  equation  UnMUnxn  =  Unf,  to  an  equivalent  equation, 
in  terms  of  solvability, 

UnRUnxn  =  Unf ,  (31) 

in  the  subspaces  [Rn]m. 


Obviously,  the  last  equation,  represents  a  discrete  system  of  equations  of  the  collo¬ 
cation  method  for  the  system  of  SIE  (25).  From  (30)  and  the  fact  that  Cn  £  [Rn]m,  we 
conclude  that  if  ( n(t)  is  the  solution  of  the  equation  (31),  then  yn(t),  defined  by 

(Pynm  =  27rl(z,  -  1)!  /  Cn(r)  {(r  ~  t)U^  log  (*  ~  l)  +  E  *kTv-k-ltk |  dr; 

(Qyn)(t)  =  27r^_  1);  J  Cn(r)r~^|(r-t)^-1log  ^1  -  0  +  '^2  PkTv~k~1tk^  dr 

(32) 

is  the  discrete  solution  of  the  system  UnMUnxn  =  Unf  and  vice  versa.  As  we  have 
mentioned  earlier,  the  vector  functions  yn(t)  are  determined  uniquely  in  terms  of  (n(t) 
according  to  formulae  (32). 

It  follows  that  if  the  equation  (31)  has  a  unique  solution  £n(i)  in  the  subspace  [Rn]m, 
then  the  equality  yn{f)  =  xn{t)  is  true.  Therefore  xn{t)  has  been  uniquely  determined. 

The  discrete  system  of  equation  (31)  is  obtained  from  the  collocation  method  ap¬ 
plied  to  (25).  We  will  show  that  for  this  system,  all  conditions  of  Theorem  8.3  from 
[23]  are  satisfied.  From  the  condition  3  of  Theorem  4  and  formula  (26),  we  have 
det[C(t)  ±D(t)]  7^  0  for  any  t  E  T,  which  coincides  with  the  condition  3  of  Theorem  8.3 
from  [23].  From  the  equality 


[C(t)  -  DitT^Cit)  +  D(t)]  =  UB-\t)Av{t) 


and  the  condition  4  of  Theorem  4,  we  conclude  the  left  partial  indexes  matrix  functions 
of  [C(t)  —  14(t)]_1[C'(t)  +  D(t )]  has  to  be  zero.  This  is  the  condition  4  of  Theorem  8.3 
in  [23]. 

From  Lemma  3  and  the  conditions  3,  4,  6  listed  in  Theorem  4,  the  reversibility 
of  operator  R  :  [Lp(T)]m  — ►  [Lp(T)]m  follows.  Note  that  the  conditions  1,  2,  5,  and 
7,  coincide  with  the  conditions  for  Theorem  8.3  from  [23].  Therefore,  if  n  is  large 
enough,  the  system  (31)  has  a  unique  solution  Q. n(t)  £  [Rn]m-  Hence,  the  linear  system 
UnMUnxn  =  Unf ,  and  the  system  of  linear  algebraic  equations  (21)  has  a  unique 
solution  if  n  is  large  enough. 


11 


From  Theorem  8.3  in  [23],  we  know  that 


IIC  -  Cn  1 1  [Lp\  m  <  O  ^  +  0(u(f\  -))  +  O  -)^  .  (33) 

From  (24)  and  (30)  we  obtain 

(PxYu\t)  =  ( PQ(t )  and  (Qx)M(t)  =  t~v  •  (QC)W- 
Therefore  we  have 

(Pxn)^(t)  =  (PCn)W, 

(Qsn)M  =  r"  (QCn)(i). 


We  proceed  to  get  an  error  estimate 

I  I  O'*  np  I  I  -  I  I  rp  (^0  /"Y»  (^0  I  I  r  t-  ^ 

llx  xn||p,*'  —  ||x  ll[Lp]m 

<  ll^(C  -  Cn)|| [Lp]rn  +  || t~UQ((  ~  Cn)\\[LPU 

<  II-PII  •  IIC- Cn||[Lp]m  +  \\t~V\\Lp  •  IIQII  •  IK  -  Cn\\[Lp]m 

<  (iiP|i  +  inMiQii)iie-enii 

=  (iiP|i  +  ciiiQii)ne-enii. 


Here  we  have  used  the  inequality 


(34) 


From  the  previous  inequality  (34)  and  with  the  help  of  (33)  we  obtain  (29).  That 
completes  the  proof  of  the  theorem. 

The  result  in  the  convergence  theorem  assume  that  the  integrations  in  the  system 
of  linear  algebraic  equations  (21)  can  be  calculated  exactly.  However,  in  practice,  those 
integrations  are  evaluated  using  the  mechanical  quadrature  rule.  We  introduce  the 
following  convergence  theorem  if  the  integrations  are  approximated  by  the  mechanical 
quadrature  rule  [23]. 


Theorem  5  Assume  that  all  the  conditions  1-7  in  Theorem  4  are  satisfied.  Then,  if 
n  is  large  enough,  the  system  of  linear  algebraic  equations  (23)  has  a  unique  solution 
(k,  k  =  —n,  n.  The  approximate  solution  of  (19)  converges  to  the  exact  solution  x(t) 
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o  M 

of  SIDE  (l)-(2)  in  the  norm  [Wp  ]m  as  n  — ►  oo.  Furthermore,  the  following  error 
estimate  is  true, 

II®  -  xn\\ |Pi„  =  5 n  +  0(ujT(h;  -))  +  O  ^log(2"+1)^  .  (35) 

n  \  na  ) 

Proof:  It  is  easy  to  verify  that  the  system  of  singular  integro-differential  equations 
(23)  is  equivalent  to  the  operator  equation 

Un\t  [Ar{t){PxV){t)  +  Br{t){Qx^){t) 

U= 0 

+  2^/  \  Un]  [ru+l~rK  (i,  r)]  (Px^)(r)dT 

r 

r 

which  is  equivalent  to  the  equation,  in  terms  of  solvability, 


Un  |c,(*)Cn(*)  +  D(t)(S(n)(t )  +  j  [rh(t,  t)\  ■  C n(T)dr  |  =  Unf,  (37) 

after  the  integral  representation  (30)  has  been  applied.  The  matrix- functions  C(t), 
D(t),  and  h(t,r )  are  determined  from  (26). 

The  equation  (37)  represents  the  system  of  equations  obtained  by  approximating 
the  integrals  in  the  system  of  SIE  (25)  with  the  mechanical  quadrature  rule. 

From  the  conditions  of  Theorem  5,  and  the  property  of  the  mechanical  quadrature 
rule  applied  to  SIE  (25),  it  is  easy  to  verify  that  (as  it  has  been  shown  in  the  proof 
process  of  Theorem  4)  the  conditions  of  Theorem  8.4  in  [23]  are  satisfied.  Therefore, 
from  Theorem  8.4  in  [23],  we  can  conclude  that  the  system  (37)  has  a  unique  solution 
Cn(t)  £  [Rn]m  if  n  is  large  enough.  Furthermore,  we  have  the  following  error  estimate, 


IK  -  0.11,1*.  <  O  (!2^i!)  +  O  (4)  +  OH/;  b) 
+  o(^(,!;i))+o(ly(ft;b). 


For  such  an  n,  the  equation  (36)  has  a  unique  solution  xn(t)  £  [Xn]m>  which  is  con¬ 
nected  with  Cn{t)  by  the  formulae  (32).  Here  we  have  used  the  fact  that  xn(t)  =  yn(t). 
Since  the  exact  solutions  x(t)  of  SIDE  (l)-(2)  and  ((t)  are  connected  by  the  formulae 
(27),  using  the  inequality  (38),  we  conclude  the  equality  (35).  Thus  we  have  completed 
the  proof  for  Theorem  5. 


Note  that  similar  results  for  the  SIDE  (l)-(2)  have  been  provided  in  Holder  spaces, 
for  example,  see  ([24]).  But  our  results  are  for  Lebesgue  spaces. 


13 


5  A  numerical  example 


We  present  the  result  of  an  example  for  m  =  1  to  show  the  performance  of  our  method. 
In  this  example,  we  take  the  exact  solution  as 


The  coefficient  matrix  functions  are  chosen  as  follows 

^M  =  AlW  =  I((  +  i-i)  (i  +  i), 

Bo(t)  =  Bl(t)  =  i(t+i-i)  (i-l), 

Kr(t,r)  =  t  +  r  +  1,  r  =  0,  1. 

T 

The  contour  T  is  an  ellipse  R cos(0)  +  irsin(/).  For  this  example,  R  =  3  and  r  =  2. 
The  right  hand  side  /(f)  then  is  determined  from  (1). 

In  Table  1  and  Figure  5,  we  show  the  result  using  the  numerical  method  developed 
in  this  paper.  The  collocation  points  are  taken  according  to  the  formula  (28). 


2  n 

Error 

8 

0.0749 

16 

0.0215 

24 

0.0012 

28 

2.8018e-004 

32 

6.4508e-005 

Table  1:  Error  between  the  exact  solution  x(z)  and  the  approximate  solution  xn(z)  at 
selected  points  z  for  different  n.  The  error  is  the  largest  error  in  the  magnitude  of  all 
selected  points. 


In  our  test,  the  non-collocation  points  have  been  obtained  from  formula 


z(j)  =  R  cos 


2n(j  —  1)  7T 

- 7 - -  +  —  )  +  r  i  sin 

k  16 


2n(j  -  1)  7 t 


k 


+  16  *  ’  i  + 


where  A;  is  a  natural  number.  We  have  observed  that,  we  should  take  enough  collocation 
points  to  guarantee  the  convergence. 


6  Conclusions 

In  this  paper,  we  have  proposed  a  numerical  scheme  for  solving  some  system  of  singular 
integro-differential  equations  for  arbitrary  smooth  contours  in  complex  plane.  The 
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(a)  (b) 


Figure  1:  Plot  of  the  solution  and  the  error  of  the  example  with  m  =  1  and  n  =  16. 
The  solution  are  evaluated  at  k  =  100  non-collocation  points,  (a)  The  solution  plot  at 
collocation  points.  The  approximate  and  the  exact  solutions  are  almost  identical  from 
the  plot,  (b)  The  error  plot  between  the  approximate  and  the  exact  solution. 


method  is  based  on  a  collocation  method  and  mechanical  quadrature  rule.  The  method 
has  been  proved  to  be  convergent  in  Lebesgue  spaces.  A  numerical  example  illustrated 
the  performance  of  the  new  method. 
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